import numpy as np
%matplotlib notebook
import matplotlib.pyplot as plt
import ipywidgets as widgets
from maux import *
hide_interactive_toolbars()

# interaktívny editor pre vyšetrovanie priebehu elementárnych funkcií
from jedit import editor
np.warnings.filterwarnings('ignore', category=np.VisibleDeprecationWarning) 

# nastavenie jazyka
from locale import setlocale, LC_ALL
from platform import uname
if uname()[0] == 'Linux':
    setlocale(LC_ALL, 'sk_SK.utf8')
else:
    setlocale(LC_ALL, 'sk_SK')
plt.rcParams["axes.formatter.use_locale"] = True

Exponenciálne funkcie

Exponenciálna funkcia je určená rovnicou $$y = a^x,$$ kde $a > 0, a \ne 1$ je pevne zvolené reálne číslo. Definičným oborom je množina reálnych čísel.

• Wikipédia
• Wikipedia

Dokumentácia:

• Exponents and logarithms

Príklad

Nakreslenie grafu exponenciálnej funkcie danej rovnicou $$y = \mathrm{e}^x.$$

#####
##### nakreslenie grafu funkcie
#####

#### vstupné údaje
def f(X): return np.e ** X # ufunc verzia funkcie
X = np.linspace(-4, 2, 6*10+1) # výber hodnôt nezávislej premennej pre zaujímavú časť grafu
Y = f(X) # odpovedajúce hodnoty závislej premennej

#### obrázok s jedným diagramom
fig, ax = plt.subplots()
fig.set_size_inches(6, 7) # veľkosť obrázka (východzia hodnota je 6x4)

### diagram
init_subplot(ax) # inicializácia diagramu: vytvorí sa pravoúhla súradnicová sústava
ax.set_title(r"Graf funkcie $y = \mathrm{e}^x$") # pomenovanie diagramu
ax.set_aspect('equal') # nastavenie rovnakej mierky pre obe osi
# ax.grid() # pravoúhla sieť

## graf funkcie
ax.plot(X, Y)

### archivácia obrázka
# fig.savefig("<meno súboru>.png")

### samotné zobrazenie
fig.show()

Príklad

Nakreslenie grafov exponenciálnych funkcií $$y = a^x$$ pre $a = \frac{1}{4}, \frac{1}{3}, \frac{1}{2}, 2, 3, 4$ do jedného obrázka.

Pokyny. Nakreslite grafy pre iný výber parametrov.

#####
##### grafy parametrického systému funkcií
#####

#### vstupné údaje
def f(X, a): return a ** X
X = np.linspace(-3, 3, 4*20+1)

#### obrázok s jedným diagramom
fig, ax = plt.subplots()
fig.set_size_inches(9, 9)

### diagram
init_subplot(ax)
ax.set_title(r"Grafy exponenciálnych funkcií $y = a^x$")
ax.set_aspect('equal')
# ax.grid()
ax.set_ylim(0, 5.5)
ax.set_xticks(range(-3, 3+1))
ax.set_yticks(range(-0, 5+1))

for a in [1/4, 1/3, 1/2, 2, 3, 4]:
    ax.plot(X, f(X, a), label=r"$a = {:.2g}$".format(a))

## legenda
ax.legend()

### archivácia obrázka
# fig.savefig("<meno súboru>.png")

### samotné zobrazenie
fig.show()

Príklad

To isté ako v predchadzajúcom príklade, len vykreslenie sa deje pomocou interaktívnych prvkov knižnice ipywidgets.

#####
##### grafy parametrického systému funkcií (interaktívna verzia)
#####

#### vstupné údaje
def f(X, a): return a ** X
X = np.linspace(-3, 3, 4*20+1)

#### obrázok s jedným diagramom
fig, ax = plt.subplots()
fig.set_size_inches(9, 9)

### diagram
init_subplot(ax)
ax.set_aspect('equal')
# ax.grid()
ax.set_ylim(0, 5.5)
ax.set_xticks(range(-3, 3+1))
ax.set_yticks(range(0, 5+1))

## graf funkcie y = aˣ
def plot_graph(a):
    ax.set_title(r"Graf exponenciálnej funkcie $y = {0}^{{x}}$".format(a))
    if ax.lines:
        ax.lines[0].set_ydata(f(X, a))
    else:
        ax.plot(X, f(X, a))

widgets.interact(plot_graph, 
                 a=widgets.SelectionSlider(options=[1/4, 1/3, 1/2, 2, 3, 4], value=2)
                 )

### samotné zobrazenie
fig.show()

Pokyny pre nasledujúce príklady

V nasledujúcich príkladoch budeme kresliť grafy a vyšetrovať priebeh exponenciálnych funkcií. Pretože tieto funkcie majú neohraničený definičný obor, budeme pri zostrojovaní grafu každej takejto funkcie vykreslovať len jej zaujímavú časť. Pri vyšetrovani priebehu týchto funkcií treba určiť:

• Nulové body funkcie.
• Extrémy funkcie.
• Intervaly, na ktorých je funkcia monotónna.
• Asymptoty grafu funkcie.

Preferujeme algebraickú metódu pri hľadaní význačných bodov. Súradnice význačných bodov, ktoré sú zrejmé z grafu, netreba explicitne uvádzať.

Doporučujeme rozdeliť riešenie do dvoch častí:

• V prvom obrázku nakreslite graf funkcie.
• V druhom obrázku vyšetrite jej priebeh.

#####
##### šablóna riešenia (nakreslenie grafu a vyšetrenie priebehu funkcie)
#####

#### vstupné údaje
# def f(X): return None # ufunc verzia funkcie
# X = np.linspace(None, None, None+1) # výber hodnôt nezávislej premennej pre zaujímavú časť grafu
# Y = f(X) # odpovedajúce hodnoty závislej premennej
# X1, X2 = X[X < None], X[X > None] # číslo None nepatrí do oboru definície
# Y1, Y2 = f(X1), f(X2) # odpovedajúce hodnoty závislej premennej

#### obrázok s jedným diagramom
# fig, ax = plt.subplots()
# fig.set_size_inches(6, 4) # veľkosť obrázka (východzia hodnota je 6x4)

### diagram
# init_subplot(ax) # inicializácia diagramu: vytvorí sa pravoúhla súradnicová sústava
# ax.set_title(r"Graf a priebeh funkcie $y = \cdots$") # pomenovanie diagramu
# ax.set_aspect('equal') # nastavenie rovnakej mierky pre obe osi
# ax.grid() # pravoúhla sieť

## x-ová os
# ax.set_xlim(None, None) # ohraničenie hodnôt pre os x
# ax.set_xticks(np.linspace(None, None, None+1)) # kótovanie x-ovej osi

## y-ová os
# ax.set_ylim(None, None) # ohraničenie hodnôt pre os y
# ax.set_yticks(None) # kótovanie y-ovej osi

## graf funkcie
# ax.plot(X, Y)
# color = ax.plot([], [])[0].get_color()
# ax.plot(X1, Y1, c=color)
# ax.plot(X2, Y2, c=color)

## nulové body funkcie
# ax.plot(None, 0, 'kx', label=r"nulový bod $?$")

## x-ové súradnice extrémov funkcie vrátane krajných bodov základného intervalu
# ps = [X[0], None, X[-1]]
# ps = [X1[0], None, X1[-1], X2[0], None, X2[-1]]

## extrémy funkcie
# ax.plot(ps[None], f(ps[None]), 'o', label=r"extrém v bode $?$")

## intervaly, na ktorých je funkcia rastúca
# color = ax.plot([], [], label=r"rastúca")[0].get_color()
# for i in None:
#     I = X[(ps[i] <= X) & (X <= ps[i+1])]
#     ax.plot(I, f(I), c=color)

## intervaly, na ktorých je funkcia klesajúca
# color = ax.plot([], [], label=r"klesajúca")[0].get_color()
# for i in None:
#     I = X[(ps[i] <= X) & (X <= ps[i+1])]
#     ax.plot(I, f(I), c=color)

## asymptota bez smernice
# Ay = np.linspace(None, None, None+1)
# Ax = None * np.ones(len(Ay))
# ax.plot(Ax, Ay, '--', lw=1, label=r"asymptota bez smernice v bode $?$")

## asymptota so smernicou
# Ax = X
# Ay = None * np.ones(len(Ax)) + None
# ax.plot(Ax, Ay, '--', lw=1, label=r"asymptota so smernicou $?x+?y+? = 0$ v bode $\pm\infty$")

## legenda
# ax.legend()
# ax.legend(loc='center') # umiestenie v strede diagramu
# ax.legend(loc=(None, None)) # ľavý dolný roh legendy má súradnice (None, None)

### archivácia obrázka
# fig.savefig("<meno súboru>.png")

### samotné zobrazenie
# fig.show()

#### editor
# editor(figure=fig, axes=ax, function=f, intervals=[X])
# editor(figure=fig, axes=ax, function=f, intervals=[X1, X2])

Úloha

Nakreslite grafy a vyšetrite priebeh týchto funkcií \begin{align*} y & = 1 + \mathrm{e}^{-x} \\ y & = \frac{\mathrm{e}^{x} - \mathrm{e}^{-x}}{2} \\ y & = \frac{\mathrm{e}^{x}+1}{\mathrm{e}^{x}-1} \\ y & = \frac{1}{3}2^{1-3x} + 2 . \end{align*}

Dôležitý príklad

Nakreslenie grafu funkcie $$y = \frac{\mathrm{e}^x-1}{x}.$$ Definičným oborom je množina pozostávajúca zo všetkých reálnych čísel rôznych od $0$.

#####
##### nakreslenie grafu funkcie
#####

def f(X): return (np.e ** X - 1) / X # ufunc verzia funkcie
X = np.linspace(-6, 3, 9*100+1) # výber hodnôt nezávislej premennej pre zaujímavú časť grafu
X1, X2 = X[X < 0], X[X > 0] # číslo 0 nepatrí do oboru definície
Y1, Y2 = f(X1), f(X2) # odpovedajúce hodnoty závislej premennej

#### obrázok s jedným diagramom
fig, ax = plt.subplots()
fig.set_size_inches(9, 7) # veľkosť obrázka (východzia hodnota je 6x4)

### diagram
init_subplot(ax) # inicializácia diagramu: vytvorí sa pravoúhla súradnicová sústava
ax.set_title(r"Graf funkcie $y = \dfrac{\mathrm{e}^x-1}{x}$", fontdict={'verticalalignment': 'bottom'}) # pomenovanie diagramu
ax.set_aspect('equal') # nastavenie rovnakej mierky pre obe osi
# ax.grid() # pravoúhla sieť
ax.set_xticks(range(-6, 3+1)) # kótovanie x-ovej osi

## graf funkcie
color = ax.plot([], [])[0].get_color()
ax.plot(X1, Y1, c=color)
ax.plot(X2, Y2, c=color)
ax.plot(0, 1, 'o', c=color, mfc='w') # číslo 0 nepatrí do oboru definície

### archivácia obrázka
# fig.savefig("<meno súboru>.png")

### samotné zobrazenie
fig.show()

#### editor
# editor(figure=fig, axes=ax, function=f, intervals=[X1, X2])